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ABSTRACT 

Aims. We aim to identify and characterize secondary planet eclipses in the original light curves of all published CoRoT planets using 
uniform detection and evaluation critetia. 

Methods. Our analysis is based on a Bayesian model selection between two competing models: one with and one without an eclipse 
signal. The search is carried out by mapping the Bayes factor in favor of the eclipse model as a function of the eclipse center time, 
after which the characterization of plausible eclipse candidates is done by estimating the posterior distributions of the eclipse model 
parameters using Markov Chain Monte Carlo. 

Results. We discover statistically significant eclipse events for two planets, CoRoT-6b and CoRoT-llb, and for one brown dwarf, 
CoRoT-15b. We also find marginally significant eclipse events passing our plausibility criteria for CoRoT-3b, 13b, 18b, and 21b. The 
previously published CoRoT-lb and CoRoT-2b eclipses are also confirmed. 
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1. Introduction 

Observing secondary eclipses of transiting extrasolar planets, 
where the planet disappears behind its host star, offers a window 
to study the planetary and orbital properties not attainable by the 
transit and radial velocity (RV) observations alone. Foremost, 
eclipses provide an indirect means to study planetary atmo- 
spheres. The two main measurables, planet's albedo and bright- 
ness temperature, are strongly coupled with the structure and dy- 
namics of the atmosphere, and can be used to educe information 
about the physical processes governing the atmosphere. Due to 
the complexity of the physics involved, this additional empir- 
ical knowledge is highly valuable when aiming to understand 
the planetary atmospheres via theoretical modeling. In addition 
to the information regarding the atmosphere, eclipses yield in- 
formation about the planet's orbit. The eclipse center times and 
durations allow us to constrain the orbit's eccentricity and argu- 
ment of periastron to a higher precision than with RV observa- 
tions alone. This is especially true for the planets orbiting faint 
or rapidly rotating stars, for which precise RV observations are 
difficult to obtain. 

A planet's emergent flux is a combination of reflected stellar 
light and the planet's thermal radiation. Both the reflective and 
thermal properties depend on the structure and composition of 
the planetary atmosphere. The reflectivity is dominated by the 
existense or the absence of high-altitude clouds, while the ther- 
mal radiation arises mainly from the heating due to the stellar 
irradiation, from tidal heating (Jackson et al. 2008| Leconte et al. 
2010) >, and from radiogenic heating. Although the reflected light 



normally dominates at visible wavelengths, the thermal radiation 
from a hot planet in a short-period orbit may make a major con- 
tribution at the red end of the visible spectrum (Lopez-Morales 
|& Seager|200"7||Fortney et al.|2008| l. 



The structure of the atmosphere and its reflective properties 
are tightly coupled, and definitive modeling of the emergent flux 
is complicated at the best. The composition of the highest op- 
tically thick atmosphere layers, and the presence or absence of 
absorbing gases or reflecting clouds, will determine the reflec- 
tivity of a planet. 

Broadly, modeling leads us to expect low geometric albedos 
(A g < 0.3) at visible wavelengths for the majority of hot Jupiters. 
Moderately hot (T e ff < 1500 K) cloud-free atmospheres are ex- 
pected to have small albedos due to the absorption by molecules 
and alkali metals ( |Burrows et al.||1997} |Marley et al.||1999| , 
while the hottest atmospheres (T e g > 1500 K) can have high- 
altitude silicate clouds leading to a high albedo (Sud arsky et al.| 
2000, 2003 ). However, strong absorption by TiO and VO gases 
can dominate over the reflectance from the silicate cloud layer, 
yielding very low albedos for the hottest atmospheres ( |Fortney| 
|et al.|[2008) . Such strong absorption has the potential to evoke 
a stratospheric temperature inversion, which would be observ- 
able as anomalously deep eclipses in the infrared. However, the 
high temperature alone does not assure the existense of a temper- 
ature inversion layer. For example, strong stellar UV radiation 
from active stars can destroy the compounds responsible for the 
high absorption, thus preventing the birth of an inversion layer 
( |Knutson et al.|2010| i. 

The observed eclipses thus far confirm a complicated pic- 
ture. Many of the secondary eclipse observations agree with a 
low albedo (e.g., |Charbonneau et al.||1999| IRowe et al.||2008]> 



but recent works by Kipping & Bakos (201 1 ) and Demory et al. 



(2011 



have also shown planets that likely have high albedos. 
Recently, [Coughlin & Lopez-Morales (2012) carried out a 
systematic search for eclipses in Kepler Q2 light curves, and note 
that "majority of the detected planet candidates emit more light 
than expected owing to thermal blackbody emission in the op- 
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tical Kepler bandpass, and present a trend of increasing excess 
emission with decreasing maximum effective planetary temper- 
ature." The authors present possible causes for the trend - in- 
cluding non-LTE or other thermal emission; fluorescent tran- 
sitions; internal energy generation; and erroneously identifying 
low-mass stars, brown dwarfs or blends as planets - but do not 
find a basis to favor any of these solutions over the others. 

The CoRoT planets are based on one of the highest-precision 
data sets available, and all of them have been thoroughly charac- 
terized in their respective discovery papers. Here we present the 
results of a homogeneous search for secondary eclipses in the 
light curves of all CoRoT planets that have been published until 
June 2012, namely CoRoT lb to 23b, including the brown dwarf 
15b, but excluding the yet unpublished planet 22b. 

CoRoT (Convection, Rotation, and planetary Transits, 
Baglin et al. 2009) is a space telescope developed jointly by 



the French space agency CNES in collaboration with ESA, 
Brazil, and seven European national space agencies. Launched 
on December 2006, it was the first space telescope dedicated to 
observing planetary transits. CoRoT offers two time cadences. 
The survey mode delivers data with a cadence of 512 s, created 
by stacking 16 exposures of 32 s. For the planet candidates iden- 
tified during an observing run, a rapid time sampling of 32 s is 
also available. Furthermore, data of bright (< R — 13.5mag ) tar- 
gets have been acquired with three-color photometry, while the 
fainter ones are observed only in a single passband. 

Our analysis is based on Bayesian statistics. The first part of 
the analysis, eclipse search, uses an approach built on Bayesian 
model selection (Gregory 2005 Ford et aT]|2007[ ). Similar ap- 
proaches have been successful when searching for planetary sig- 
nals from RV data (|Gregory|2007a|b"{ |Tuomi & Kotiranta|2009| 



Greg ory & Fischer|2010| l, reanalyzing planetary systems having 
multiple degenerate solutions ([Gregory 2005 ), and disproving 
proposed planetary RV signals with insufficient evidence (|Tuomi| 



201 1 Gregory |2011| l. The second part, eclipse characterization, 



uses Bayesian parameter estimation to derive the posterior dis- 
tributions for the parameters describing the planet and its orbit 
( |Ford|2005 1 |Gregory|2005) l 



Two of the CoRoT planets, C oRoT-lb ( |Barge et al.1[2008] i 
and CoRoT-2b ( Alonso et al. 2008 ), have been reported with sec- 



ondary eclipses in the CoRoT light curves ( Snellen et al. 



2009 



Earfh- 



|Alonso et al.|2009a|b[ ), that have later been confirmed by 
based observations in the IR (|Rogers et aT||2009| Gillon et al. 
lQ09j |Alonso et al.|2010| |Zhao et al.|2012) . These two planets 



are used to verify our method. 

A brief introduction to the Bayesian model selection and 
parameter estimation approaches follows in Sec. [2] . The data 
used in the analysis are summarized in Sec. [3] and the practical 
side of our analysis methods in Sec. [4] We present our results in 
Sec. [5] and discuss the detected eclipse candidates individually 
in Sec. [6] Finally, we review the main implications of our results 
and the use of Bayesian model selection in eclipse searches in 
Sec. [3 



2. Theory 

2.1. Bayesian statistics 

We start with a short overview on the Bayesian statistics, for 
an in-depth treatise of the Bayesian approach, see, for example, 
|Gregory| ( |2005l l, |Kass & Raftery| ( [T995] l, and |Robert| ( |2007] l. 

Bayes' theorem states that the posterior probability distribu- 
tion P(A\D) for a proposition A with given data D can be derived 
from the prior probability distribution P(A) and the likelihood of 



the data for the given proposition, P(D\A), as 

P(A)P(D\A) 



P(A\D) 



J P(A)P(D\A)dA ' 



(1) 



For a model selection problem, different propositions correspond 
to different models, while for parameter estimation the proposi- 
tions correspond to different values of a model-specific parame- 
ter vector 0. 

If we assume the observed light curve fluxes Fq to follow 
from a model M(0, t), where t is the exposure center time, and 
further assume independent and identically distributed (i.i.d.) er- 
rors (deviations from the model) that follow a normal distribu- 
tion with a standard deviation cr, the probability P(e,|0) for an 
error e-, = Fo(tj) - Fm(0, f,) of the light curve point i is 



P(eM = 



1 



•V2^ 



exp 



± 
2a- 1 



(2) 



and the likelihood (joint probability) for N i.i.d. errors is ob- 
tained using the multiplication rule of independent probabilities 



P(D\0) = p| 

i=l 



1 ( e] 



= (cr ^hi) 



exp 



A e 2 ^ 
Zj 2ct 2 



(3) 



Since the number of points in the CoRoT light curves is 
large, it is necessary to carry out the calculations using the natu- 
ral logarithm of the likelihood, which is now 



N , 
lnP(D\0) = --\n2n-Nlncr-x ■ 



(4) 



Finally, the logarithm of the posterior probability for a model M 
with a parameter vector 6 given the observed flux Fq is 

\nP M {0\F Q ) = lnP(0) - | ln27r- N In cr - X 2 ■ (5) 



2.2. Parameter estimation 



Parameter estimation in its simplest form boils down to finding 
the posterior maximum for the model-specific parameter vector. 
If the posterior is likelihood dominated, this equals to maximiz- 
ing the likelihood function, and, if the error is treated as constant, 
minimizing x 2 - 

For a more robust analysis, we usually want to compute the 
full posterior distributions for the parameters, from which we 
can derive different summary statistics. In the current study we 
compute the final posterior distributions using Markov Chain 
Monte Carlo (MCMC) sampling. 



2.3. Model selection 

In order to compare the ability of our competing models to ex- 
plain the data, we need to calculate the global likelihood for a 
model by marginalizing (integrating) the posterior probability 
over the model parameters as 



P(M) 



! 



P(6)P(D\0)d0. 



(6) 
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Table 1. Rough guidelines to the interpretation of the Bayes fac- 



modified 


version of a 


similar tabulation by | Jeffrey s| ((1961 


p. 




21n(B 10 ) 


610 


Evidence against Mo 


0to2 


1 to 3 


Not worth more than a bare mention 


2 to 6 


3 to 20 


Positive 


6 to 10 


20 to 150 


Strong 


> 10 


> 150 


Very strong 



Given the global likelihoods computed for each model, we can 
calculate the Bayes factor (Kass & Raftery 1995, Gregory [ 2005 



Jeffreys|1961| l of a model M, over model Mj as 



B,j = P(Md/P(Mj), 



(7) 



and given the Bayes factors, we can calculate the normalized 
model probabilities for all the models as 



P(Mi) = 



Bio 



2 j=0 Bjo 

With two models, the probability for the model M\ is given by 

1 



(8) 



P(M t ) 



1 + 1/B W ' 



(9) 



and for the model M we get f(M ) = 1 - P(M\). 

The Bayes factors can also used directly to give "a rough 
descriptive statement about standards of evidence" (Kass & 
|Raftery|[T9 95), and we reproduce the guidelines for the direct 
interperation as given by |Kass & Raftery | in Table [T] 

2.4. The role of prior distributions 

The prior distributions (priors) encode our knowledge about the 
plausible parameter values and relative model probabilities. In 
parameter estimation, the importance of a prior depends on the 
relative amounts of parameter-specific information encoded by 
the prior and information incorporated into the data. If we al- 
ready have an estimate for a parameter, e.g., from previous re- 
search, we can use the estimate to construct an informative prior 
that may have a substantial role in the posterior distribution. In 
the absence of a prior estimate, we use an uninformative prior 
that aims to maximize our uncertainty regarding the parameter 
value, and the posterior will be dominated by the likelihood. 

We use two uninformative priors in our work. The Uniform 
prior, 

P(x; a, b) = 1 lib - a), a < x < b, (10) 
is used for location parameters, and the Jeffreys' prior, 

1 



P(x; a, b) 



xln{b/a) 



< a < x < b, 



(11) 



is used for scale parameters (for a more detailed explanation of 
location and scale parameters, see, e.g., [Gregory 2005). The pa- 
rameters a and b define the span of the prior. For the MCMC 
characterization we also use a Normal prior 

"^--^s^Hc^f)- <i2> 

where fi is the mean and <x is the standard deviation of the distri- 
bution. 



As long as the prior distribution encompasses the likelihood 
function, the exact form of an uninformative prior has often only 
a small effect in parameter estimation. However, in model selec- 
tion, the prior widths reflecting the state of our ignorance can 
affect the global likelihood significantly. 

The greater our uncertainty about a parameter value is, the 
wider the prior will be. Since priors are normalized probability 
distributions, a wider prior yields a smaller posterior point value, 
and a smaller marginal likelihood. 

Thus, if our prior distributions are wider than the likelihood 
distributions, each additional model parameter adds a penaliza- 
tion factor to the global likelihood. These factors are multiplica- 
tive, and work to penalize complex models against simpler ones. 
For a complex model to be considered better than a simpler one, 
the additional complexity must be supported by a sufficiently 
improved likelihood. 



3. Data 

For the analysis of CoRoT-lb to 21b and 23b, we use the lat- 
est versions of the CoRoT N2 light curves available publicly 
from the IAS Data CentetQas of June 2012. We use the chro- 
matic light curves when available, but include from them only 
the red channel for the analysis to reduce the number of jumps 
from cosmic ray hits. We use both the 32 and 512 second time 
cadences (when available) by assigning the two cadences a sep- 
arate mean point-to-point scatter (error) estimate in the analysis, 
and exclude the phase-span near the primary transit completely 
from the analysis. 

The stellar and planetary properties used in constructing the 
priors are gathered from the latest publications discussing the 
systems. The adopted planet properties are listed in Tableland 
the host star properties in Table[3j both sorted by the presence of 
eclipse events, as indicated in Sec. [6] 



4. Analysis 

4.1. Overview 

Our analysis consists of three main steps. First, we identify the 
eclipse candidates and assess their significances. Next, we carry 
out tests for the plausible candidates to verify that they are not 
produced by singular events (jumps) or correlated noise. Finally, 
we characterize the plausible eclipse candidates that pass our 
tests by calculating the marginal posterior probabilities for the 
eclipse model parameters. 

Our code combines Fortran2003 for the numerically inten- 
sive computations and Python^] for the high-level functionality. 
We use simple Monte Carlo (MC) importance sampling to ob- 
tain the global likelihood estimates, and Markov Chain Monte 
Carlo (MCMC) methods to obtain the parameter posterior dis- 
tributions. Parallelization is carried out using OpenMP and MPI. 
We use the NumP>P1 SciPy 3 , MPI4P>P1 matplotlitf| PyFITS^ 
and emcee (Foreman-Mac ke y et al.|2012| l -libraries for the data 
analysis, MPI parallelization, and data IO. 

We use the Mandel & Agol (2002) transit shape model with- 
out limb darkening to model the eclipse. The model is imple- 
mented in Fortran2003 and uses OpenMP for parallelization 
within a single computing node. 



1 


idoc- corot . ias . u- psud . f r| 2 


www . python . org 


3 


www.scipy.org, 4 mpi4py.scipy.org, 5 


matplotlib.org 


6 


www. stsci .edu/institute/software_hardware/pyfits 
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Table 2. Principal parameters of the analyzed CoRoT planets 
sorted as a function of the significance of the found eclipse 
events. Listed are the orbit period p , scaled semi-major axis # a 
(a/R+), planet mass Mp, and planet radius Rp. 



Planet 



0„ [d] 



M P [M Jup 



Ref. 



Systems with significant eclipse events 



CoRoT-lb 

CoRoT-2b 

CoRoT-6b 

CoRoT-llb 

CoRoT-15b 



CoRoT-3b 
CoRoT-13b 
CoRoT-18b 
CoRoT-21b 



CoRoT-4b 

CoRoT-5b 

CoRoT-7b 

CoRoT-8b 

CoRoT-9b 

CoRoT-lOb 

CoRoT-12b 

CoRoT-14b 

CoRoT-16b 

CoRoT-17b 

CoRoT-19b 

CoRoT-20b 

CoRoT-23b 



1.51 
1.74 
8.89 
2.99 
3.06 



4.9 
6.7 
17.9 
6.9 
6.7 



1.03: 
3.31 : 
2.96: 
2.33: 

63.30: 



0.12 
0.16 
0.34 
0.34 
4.10 



1.49 ±0.08 
1.47 ±0.03 
1.17 ±0.04 
1.43 ±0.03 

i 19+0.30 
1 - lz '-0.15 



Systems with marginal eclipse events 
4.26 7.8 21.77 ±1.00 1.01 ±0.07 
4.04 10.8 1.31 ±0.07 0.89 ±0.01 
1.90 6.4 3.47 ±0.38 1.31 ±0.18 
2.72 4.6 2.26 ±0.31 1.30 ±0.14 



1 

2,3 
5 
6 



4 
7 
9 
10 



Systems 

9.20 

4.04 

0.85 

6.21 
95.27 
13.24 

2.83 

1.51 

5.35 

3.77 

3.90 

9.24 

3.63 



without identified eclipse 



17.4 
9.0 
4.3 
17.6 
93.0 
31.3 

4.9 
11.2 

6.2 

6.7 
19.0 

6.9 



0.72 ± 0.08 
47+ 05 

"•^ ' -0.02 

0.02 
0.22 ± 0.03 
0.84 
2.75 
0.92 
7.60 
0.54 
2.43 
1.11 
4.24 
2.80 



1 q+o.06 

• -0.05 



0.07 
0.16 
0.07 
0.60 
0.09 
0.30 
0.06 
0.23 
0.30 



events 
1. 

1.33 ±0.05 

0.15 
0.57 ± 0.02 
0.94 ± 0.04 
0.97 ± 0.07 
1.44 ±0.13 
1.09 ±0.07 
1.17 ± 0.15 
1.02 ±0.07 
1.29 ±0.03 
0.84 ± 0.04 
1.05 ±0.13 



11, 



14, 



12 
13 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 



References. (1) Barge et al. ( 2008 ); (2) Alonso et al. ( 2008 1; (3) Bouchy 
|et al.| <|2008>; (4) |Deleuil et al.| <[2008|>; (5 ) JFridjund et al.H2010>; 
(6)|Gandolfi et al-H2010|; (7)|Cabrera et al.H2010>; (8)IBouchy et al. 
(|2010| > ; (9)|Hebrard et al.fpOl 1); (10)|Patzold et al.|(2012|>; (1 l)|Aigrain 
[eTldTp008ir; (12) IMoutou et al.| <2008|>; ( 13) |Rauer et al.|T2 009 ); 
(14)|L6ger etaT] (|2009|>; (15)|Bruntt etal] (|2010b; (16)|Borde eTal 
l|2010); (17)|Deeg et al.|(|20T0); (18)|Bonomo et al.|p010); (19 ) |Giiion 
|et al.| p010|i; (20) |Tingley et al.| j20Tlb|>; (21)|011ivier et al.||2012>; 

~p011|>; (23)[Guenther et al.| l |2012| >; (24) |Deleuil| 



(22) Csizmadia et al. 



|et al.| < 12012) ; (25) |Rouan eTaLl<2012> . 



The computations were carried out using the Diodo clus- 
ter located at the Instituto de Astroffsica de Canarias and the 
LaPalma supercomputer located at the La Palma Astrophysics 
Centre. 



4.2. Models and Parameterization 

We consider two models: with and without an eclipse (Mi and 
Mo, respectively). Table [4] offers an overview of the model pa- 
rameters and their priors. With Mo being a submodel of Mi, the 
models share parameters, which allows us to sample both of the 
model likelihood spaces simultaneously. 

The no-eclipse model, Mo, attempts to explain the data with 
a constant zeropoint value, Ofo, and mean point-to-point scatter 
estimates for the light curves of both time cadences (Oo-i and 6^, 
respectively). Thus, Mo has 2-3 free parameters, depending on 
whether both time samplings are available for a given CoRoT 
target. 

The eclipse model, Mi , includes the parameters of Mo, and 
introduces eight additional parameters. However, only three of 



these additional parameters have posterior distributions domi- 
nated by the likelihood from the data. The five other additional 
parameters have prior-dominated posteriors, and are included to 
propagate the uncertainties in their estimates into the analysis. 

Of the prior-dominated parameters, the radius ratio, period, 
scaled semi-major axis, and inclination have priors based on 
published values that were estimated from the primary transits. 
The contamination factor gives the amount of third light inside 
the CoRoT aperture mask ( |Deeg et al.|2009"l|Deleuil et al.|2009) , 
and has a prior based on the published contamination values. 

Of the likelihood-dominated parameters, the planet-star flux 
ratio Of gives the ratio of the surface brightness of the planet 
to the surface brightness of the star, and is related to the 
eclipse depth AF and planet-to-star surface area ratio 8^ as 
8f = The eccentricity and argument of periastron de- 

termine the eclipse center and duration. Their posterior proba- 
bilities are constrained by radial velocity orbservations, but the 
constraints are often wide. 

Both models also include two derived parameters: the eclipse 
center and duration, T e and D e , respectively. The dependency 
of Mo on these two parameters is due to our chosen detrending 



approach, and is explained in more detail in Sec. 4.4 



4.3. Parameter priors 

As mentioned previously, priors have an important role in the 
Bayesian model selection. Since the eclipse model, M\, has eight 
parameters more than the no-eclipse model, Mo, it would be ex- 
pected that the penalization against M\ would be strong. 

However, the five prior-dominated M\ parameters are closer 
to fixed constants than to free parameters, and the penalization 
by most of them has a relatively small effect (the likelihood 
within the prior limits is close to a constant). For these parame- 
ters, the additional information in the eclipse is very small com- 
pared to the information from the primary transit and radial ve- 
locity observations. The main motivation for including the prior- 
dominated parameters comes from the added robustness. The pa- 
rameters are included as nuisance parameters to be marginalized 
over in order to propagate their uncertainties to the global likeli- 
hoods. 

Excluding the unknown flux ratio 8f, all the Mi parameters 
have uniform priors centered around the latest published values 
and half-widths corresponding to their l<x uncertainties. In many 
cases the eccentricities are poorly constrained by the radial ve- 
locity measurements, and we set the priors to span the eccentric- 
ity space from zero to the published maximum value. 

The zeropoint and point-to-point scatter are handled some- 
what differently than the other parameters. The posteriors of 
these three parameters were found to be close to normal when 
the other parameters are held fixed. This allows us to use 



Laplace's approximation (see, for example, Gregory 2005 1 to 
obtain an analytical likelihood estimate marginalized over the 
three parameters for each MC sample, and reduce the sampling 
variance without sacrificing accuracy. 

4.4. Light curve detrending 

Since the depth of an eclipse is on the order of percents of 
the depth of a transit, great care must be taken to remove the 
signals not related to an eclipse. CoRoT light curves contain 
four major sources of correlated noise: stellar variability on long 
time-scales, short time-scale variability from stellar granulation 
( |Aigrain et al.|[2009] >, random jumps due to cosmic ray hits 
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Table 3. Basic host star properties for the analyzed CoRoT planets. 



Star Spectral type V T t [K] M+ [M Q ] R+ [R ] logg [Fe/H] Contamination [%] Ref. 



Systems with significant eclipse events 



CoRoT-1 


GOV 


13.6 


5950 ± 


150 


0.95 ± 


0.15 


1.11 + 0.05 


4.25 ± 


0.30 


-0.30 


±0.25 


~ 


1 


CoRoT-2 


G7V 


12.6 


5625 ± 


120 


0.97 + 


0.06 


0.90 ± 0.02 


4.30 + 


0.20 


0.00 


±0.10 


5.6 


2,3 


CoRoT-6 


F9V 


13.9 


6090 ± 


50 


1.05 ± 


0.05 


1.03 + 0.03 


4.44 + 


0.23 


-0.20 


±0.10 


2.8 


5 


CoRoT-11 


F6V 


12.9 


6440 + 


120 


1.27 ± 


0.05 


1.37 + 0.03 


4.22 + 


0.23 


-0.03 


±0.08 


13.0 


6 


CoRoT-15 


F7V 




6350 ± 


200 


1.32 ± 


0.12 


1 46 +031 


4.30 + 


0.20 


0.10 


±0.20 


1.9 


8 



Systems with marginal eclipse events 



CoRoT-3 


F3V 


13.3 


6740 ± 


140 


1.37 ± 


0.09 


1.56 ± 


0.09 


4.22 ± 


0.07 


-0.02 


±0.06 


8.2 


4 


CoRoT-13 


GOV 


15.0 


5945 ± 


90 


1.09 ± 


0.02 


1.01 ± 


0.03 


4.30 ± 


0.10 


0.01 


±0.07 


11.0 


7 


CoRoT-18 


G9V 


15.0 


5440 ± 


100 


0.95 ± 


0.15 


1.00 ± 


0.13 


4.40 ± 


0.10 


-0.10 


±0.10 


2.0 


9 


CoRoT-21 


F8IV 


16.1 


6200 ± 


100 


1.29 ± 


0.09 


1.95 ± 


0.21 


3.70 ± 


0.10 


0.00 


±0.10 


8.5 


10 



Systems without identified eclipse events 



CoRoT-4 


F8V 


13.7 


6190 ± 


60 


1.16 ± 


0.03 


1.17 ± 


0.03 


4.41 


±0.05 


0.05 


±0.07 


0.3 


11, 12 


CoRoT-5 


F9V 


14.0 


6100 ± 


65 


1.00 ± 


0.02 


1.19 ± 


0.04 


4.19 


±0.03 


-0.25 


±0.06 


8.4 


13 


CoRoT-7 


K0V 


11.7 


5250 ± 


60 


0.91 ± 


0.03 


0.82 ± 


0.04 


4.57 


±0.05 


0.12 


±0.06 


0.5 


14, 15 


CoRoT-8 


K1V 


14.8 


5080 ± 


80 


0.88 ± 


0.04 


0.77 ± 


0.02 


4.58 


±0.08 


0.30 


±0.10 


0.7 


16 


CoRoT-9 


G3V 


13.7 


5625 ± 


80 


0.99 ± 


0.04 


0.94 ± 


0.04 


4.54 


±0.09 


-0.01 


±0.06 


2.5 


17 


CoRoT-10 


K1V 


15.2 


5075 ± 


75 


0.89 ± 


0.05 


0.79 ± 


0.05 


4.65 


±0.10 


0.26 


±0.07 


5.5 


18 


CoRoT-12 


G2V 


15.5 


5675 ± 


80 


1.08 ± 


0.08 


1.10 ± 


0.10 


4.52 


±0.08 


0.16 


±0.10 


3.3 


19 


CoRoT-14 


F9V 


16.0 


6035 ± 


100 


1.13 ± 


0.09 


1.21 ± 


0.08 


4.35 


±0.15 


0.05 


±0.15 


7.0 


20 


CoRoT-16 


G5V 


15.6 


5650 ± 


100 


1.10 ± 


0.08 


1.19 ± 


0.14 


4.36 


±0.10 


0.19 


±0.06 


2.3 


21 


CoRoT-17 


G2V 


15.5 


5740 ± 


80 


1.04 ± 


0.10 


1.59 ± 


0.07 


4.40 


±0.10 


0.00 


±0.10 


8.0 


22 


CoRoT-19 


F9V 


14.8 


6090 ± 


70 


1.21 ± 


0.05 


1.65 ± 


0.04 


4.07 


±0.03 


-0.02 


±0.10 


0.3 


23 


CoRoT-20 


G2V 


14.7 


5880 ± 


90 


1.14 ± 


0.08 


1.02 ± 


0.05 


4.20 


±0.15 


0.14 


±0.12 


< 0.6 


24 


CoRoT-23 


GOV 


15.6 


5900 ± 


100 


1.14 ± 


0.08 


1.61 ± 


0.18 


4.30 


±0.20 


0.05 


±0.10 


7.2 


25 



Refe rences. (l)|Barge et aLlpOOg]); (2)|Alonso eTaIlp008l; (3)|Bouchy et al.|p008l; (4)|Deleuil et al.| p008>; (5)|FndTund et al.| (f20Tol> 
(6)|Gandolfi et al.H2010|>; (7)|Cabrera et al.|l|2010); (8)|Bouchy et aT]l|2010|>; (9)|Hebrard et al.lpOl 1); (10)|Patzold et al.H2012[>; (1 l)|Aigrain et al. 



(12)|Moutou et al.|p008> ; (13)|RauerltaT1(|2009|); (14)|Leger eTaLlj2609>; (15)|Bruntt etaL|j20iO|>; (16)|Borde et~at1 ( |20lb|>; (17)|Deeg 



|et al.|d2010[ > ; (18)|Bonomo et al.|l|20T0); (19)|Gillon et al.| <|2010| >; (20) Tingley et alT|l |2011b| >; (21) |011ivier et al.| < |2012| >r(22) Csizmadi a et al. 
POII) ; (23) [Guenther et al.|P012| >; (24) |Deleuil et al.| ( |2012) ; (25) |Rouan et al.| l |2012| >. 



(Pinheiro da Silva et al.||2008), and a quasi-periodic signal re- Aigrain et al.|[2009). The signal is quasi-periodic over the pe 



lated to the satellite's orbit ( jPinheiro da Silva et al . 2008 ; Aigrain 
|et al.|2009) . 

We chose a detrending approach in which the detrending is 
carried out separately for each planetary period, and depends on 
the period, eclipse center and eclipse duration. Thus, we need to 
detrend the data separately for each MC integration sample. This 
increases the computational burden (the detrending is one of the 
most computationally intensive tasks in the analysis), but allows 
us to average over the possible effects caused by the detrending. 

The primary source of light curve variability is usually the 
star itself. Fortunately, since the time scales for stellar variability 
are significantly longer than the eclipse duration, this variability 
is also the easiest external signal to remove: the trends can be 
modeled locally with a polynomial fitted to the out-of-eclipse 
(OE) data around each eclipse. 

The second major noise signal is due to high-energy parti- 



cles impacting the detector ( Pinheiro da Silva et al.| |2008). The 
impacts cause sudden jumps in the flux level, followed by either 
gradual or abrupt return close to the original level. While meth- 



ods have been developed to correct for these jumps (|Mislis etaL 



2010), we decide to simply exclude the orbital periods with iden- 
tifiable jumps within the eclipse search range. This is a prefer- 
able approach since even a small error in the correction could 
still yield a signal strong enough to mimic an eclipse. 

The third noise signal comes from the crossing of the South 
Atlantic Anomaly (SAA) and the satellite's entering and ex- 
iting from the Earth's shadow (Pinhe iro da Silva et al.||2008| 



riod of the satellite's orbit, has relatively small amplitude, and 
does not change significantly over the duration of an eclipse. We 
first reduce the signal by excluding the points flagged by the N2 
pipeline as having been obtained during the SAA crossing. The 
residual signal is modeled and removed using a periodic spline 
fitted to the OE flux around each eclipse folded over the satel- 
lite's period. 

Finally, while the variability due to stellar granulation has a 
low amplitude, it can still be orders of magnitude greater than an 
eclipse signal (especially, K-dwarfs and giants have been shown 
to feature variations on a 0.5 mmag level on 2 h time scales, 
Aigr ain et al.| [2009 ). Since the variability is on the same time- 
scale as the eclipses, it cannot be removed safely. The Bayes 
factor mapping of light curves with substantial short time-scale 
variation usually shows several high-probability peaks with un- 
realistic eclipse depths. These peaks from single correlated-noise 
events can often be identified at the later stages of the analysis, 
as explained in Sec. |4.6| but may also render the analysis of faint 
targets infeasible. 

4.5. Eclipse search 

The eclipse search is based on Bayesian model selection. We 
carry out the most numerically intensive task, the likelihood 
sampling, separately from the rest of the analysis. After we have 
the likelihood samples, the computation of posterior probabil- 
ities, marginal posteriors, global likelihoods and Bayes factors 
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Table 4. Model parameterizations. Total number of parameters depends on whether the light curves include both time cadences. 
The third and fourth columns list whether the parameter is included in Mq and M\ models, respectively. The prior column lists the 
type of the prior used for the parameter. Marginalization lists the method used to compute the likelihood, MC stands for Monte 
Carlo sampling and LA for Laplace approximation. The effects from the zeropoint correction are small, but are still included for 
thoroughness. 





Notation 


M 


Mi 


Prior 


Marginalization 


Notes 


Planet-star flux ratio 


Oi 




X 


Jeffreys' 


MC 


Likelihood dominated 


Planet-star radius ratio 


0, 




X 


Uniform 


MC 


Prior dominated 


Period 


% 




X 


Uniform 


MC 


Prior dominated 


Scaled semi-major axis 


e a 




X 


Uniform 


MC 


Prior dominated 


Eccentricity 


e e 




X 


Uniform 


MC 


Likelihood dominated 


Argument of periastron 


o a 




X 


Uniform 


MC 


Likelihood dominated 


Inclination 






X 


Uniform 


MC 


Prior dominated 


Contamination" 


o c 




X 


Uniform 


MC 


Prior dominated 


Zeropoint 


Ota 


X 


X 


Uniform 


LA 




Long cadence error 


6<rl 


X 


X 


Jeffreys' 


LA 




Short cadence error 


Oas 


X 


X 


Jeffreys' 


LA 




Eclipse duration 


£>e 


X 


X 


Uniform 


MC 


Derived from 6> c and G bJ 


Eclipse center time* 




X 


X 


Uniform 


MC 


Derived from 6> c and U 



Notes. <a) Contamination measures the amount of third light inside the CoRoT aperture mask. (fc) Eclipse center time is the time relative to 
equidistance between primary transits. 



is fast. Also, this separation allows us to evaluate the effect that 
different priors have on the posteriors and on the Bayes factors. 

We start by calculating a set of likelihood samples for the 
two competing models, Mq and M\, using Monte Carlo (MC) 
importance sampling combined with the Laplace approximation 
(LA). The MC sampling ranges are based on the published pa- 
rameter values with their errors. The Laplace approximation for 
the likelihood marginalized over the zeropoint and point-to-point 
scatter is calculated for each MC sample. The approximation is 
calculated by first maximizing the likelihood as a function of the 
LA parameters, then computing the covariance matrix in these 
dimensions around the maximum, and finally using the obtained 
covariances to calculate the integral (for a detailed example, see 
Gregory|2005l p. 291). 

Given the likelihood samples for both models, we calculate 
the Bayes factors for a set of eclipse centers, C. The set spans the 
whole ecplice center space (determined by the 9 p , 9 e , and M pri- 
ors) linearly. We compute the global likelihoods with a constant- 
width uniform prior on the eclipse center time. For each r e € C, 
the prior is centered around T e , and has a width of one tenth 
of the transit duration. The method generates a one-dimensional 
Bayes factor map, from which the eclipse candidates can be 
identified (for example, see Fig. [T}. The approach corresponds 
to making tic model comparisons, where tic is the number of 
suggested eclipse centers, where each comparison differs by the 
priors imposed on the parameters. 

4.6. Plausibility tests 

After the initial identification of the eclipse candidates, further 
steps must be taken to assess their significance. The most impor- 
tant criterion we use for testing the plausibility of a candidate 
is the "splitting test", where we require that the signal must be 
found from separate subsections of the light curve. A jump by 
a high-energy particle hit or correlated noise , e.g., from stellar 
activity, can both generate dips that mimic an eclipse signal, but 
these signals arise from single events. 



To separate signals arising from single events from plausible 
eclipse signals, we map the Bayes factors for each CoRoT target 
featuring a significant signal for three light curve subsets. We 
remove one thirds of the total light curve from each subset, and 
require that the eclipse signal is visible for all the subsets. This 
approach allows us to reject most of the false eclipse candidates 
arising from single events, but cannot be used to test the weakest 
signals. 

4.7. Eclipse characterization 

After identifying the most plausible eclipse candidates, we carry 
out MCMC simulations to obtain the posterior distributions for 
the Mi parameters. We originally used our own implementation 
of the Metropolis-Hastings sampler, but changed to use emcee 
by [Foreman-Mackey et al. ( 2012| ), a freely available Python im- 
plementation of the Affine Invariant Markov Chain sampler with 
excellent convergence properties (Goodm an & Weare||2010| . 
The change of the sampler did not alter the results, but reduced 
the computation time significantly due to the improved poste- 
rior sampling properties (smaller number of MCMC steps were 
needed in order to obtain independent samples). 

We change our parameterization slightly for the character- 
ization step. First, we assign proper priors for the zeropoints 
and errors (uniform and Jeffreys', respectively). Second, since 
we already have an estimate for the eclipse center from the 
Bayes factor mapping, r max(Bl0 ), we assign a normal prior N(fi = 
7maxOB 10 )> °~ = 0.2 h) on the eclipse center. We also add a prior 
on the primary transit duration, which allows us to constrain the 
(6> e , 6> M )-space. 



5. Results 

We summarize our results for the planets with eclipse candidates 
in Tableland Figs.[T]to[5] The first two columns in Table[5]list 
the Bayes factors and the corresponding Mi posterior probabili- 
ties from the model comparison analysis. The next four columns 
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Table 5. Main results form the combined Bayesian analysis for the planets with detected eclipse candidates. Listed are the maximum 
Bayes factor B 10 , and the Mi posterior probability P(M\ ); MCMC estimates for the eclipse center phase, orbit eccentricity, flux ratio, 
and eclipse depth; and estimates for the brightness and equilibrium temperatures. 



Planet 


Bio 




P(M0 [%] 


Phase Eccentricity Flux ratio [%] 


Depth [% ] 




[K] 


T 


[K] 












Previously reported planets 












CoRoT-lb 


7.6 ± 


0.2 


88.4 ± 


0.3 


0.495 ± 0.002 0.04 ± 0.03 0.99 ± 0.39 


0.20 ± 0.08 


1580 


-2500 


1830- 


-2430 


CoRoT-2b 


2.0 ± 


0.0 


66.1 ± 


0.2 


0.501 ±0.002 0.03 ±0.03 0.25 ±0.16 


0.07 ± 0.05 


1430 


-2110 


1460- 


-1930 












Statistically significant new eclipse candidates 












CoRoT-6b 


10.6: 


±0.7 


91.4 ± 


0.4 


0.533 ±0.000 0.06 ±0.01 1.96 ±0.69 


0.27 ±0.10 


2230 


-2840 


970- 


1270 


CoRoT-llb 


(5.5 = 


t 1.1) x 10 4 


100.0 : 


±0.0 


0.558 ± 0.002 0.35 ± 0.03 3.15 ± 0.58 


0.36 ± 0.07 


2580 


-3020 


1650 


-2180 


CoRoT-15b 


27.2: 


±2.9 


96.5 + 


0.4 


0.495 ± 0.001 0.08 ± 0.05 22.00 ± 6.50 


1.37 ±0.40 


3470 


-4710 


1580- 


-2340 












Statistically marginal new eclipse candidates 












CoRoT-3b 


2.3 ± 


0.0 


69.7 ± 


0.3 


0.509 ±0.001 0.06 ±0.06 1.85 ±1.15 


0.08 ± 0.05 


1940 


-3070 


1610- 


-2170 


CoRoT-13b 


3.3 ± 


0.1 


77.0 + 


0.5 


0.483 ± 0.001 0.08 ± 0.04 3.04 ± 1.47 


0.25 ±0.12 


2150 


-3080 


1220- 


-1610 


CoRoT-18b 


3.1 ± 


0.1 


75.9 ± 


0.4 


0.469 ± 0.002 0. 10 ± 0.04 3.93 ± 1 .67 


0.71 ±0.30 


2130 


-3020 


1440- 


- 1950 


CoRoT-21b 


1.7 ± 


0.0 


63.1 ± 


0.3 


0.474 ± 0.002 0.05 ± 0.02 5.99 ± 3.29 


0.27 ±0.15 


2210 


-3600 


1930- 


-2600 



Table 6. Eclipse ephemerides for the detected eclipse candi- 
dates. The periods are taken from the references listed in Table 




Planet 


Eclipse center 


Period 




[HJD-2454000] 


[d] 


CoRoT-lb 


160.2001 ± 0.0030 


1.50896 


CoRoT-2b 


238.4089 ± 0.0035 


1.74210 


CoRoT-3b 


285.3050 ± 0.0043 


4.25680 


CoRoT-6b 


600.3510 ± 0.0009 


8.88659 


CoRoT-llb 


599.3498 ± 0.0060 


2.99433 


CoRoT-13b 


792.7581 ± 0.0040 


4.03519 


CoRoT-15b 


755.0757 ± 0.0030 


3.06036 


CoRoT-18b 


1322.6153 ± 0.0038 


1.90007 


CoRoT-21b 


400.3197 ± 0.0054 


2.72474 


Table 7. Expected flux ratios, eclipse depths, 


and equilibrium 


temperatures for the planets without detected eclipse candidates. 


Shown are 


the 95% credible intervals for the expected Of, AF, 


and r eq when the planet is assumed to radiate as a blackbody. 


Planet 


Flux ratio [%] Depth [%o] 


T eq [K] 


CoRoT-4b 


0.01-0.15 0.00-0.02 


1010-1350 


CoRoT-5b 


0.05-0.58 0.01 -0.11 


1380- 1850 


CoRoT-7b 


0.44-3.13 0.00-0.01 


1720-2350 


CoRoT-8b 


0.00-0.14 0.00-0.01 


820- 1110 


CoRoT-9b 


0.00-0.01 0.00-0.00 


390- 530 


CoRoT-lOb 


0.00-0.05 0.00-0.01 


610- 840 


CoRoT-12b 


0.07-0.78 0.01-0.18 


1380- 1860 


CoRoT-14b 


0.46 - 2.92 0.04 - 0.25 


1850-2560 


CoRoT-16b 


0.02 - 0.43 0.00 - 0.04 


1110-1610 


CoRoT-17b 


0.15- 1.33 0.01 -0.06 


1550-2110 


CoRoT-19b 


0.15-1.13 0.01-0.07 


1590-2140 


CoRoT-20b 


0.00-0.13 0.00-0.01 


910-1240 


CoRoT-23b 


0.11- 1.27 0.01-0.06 


1500-2130 



show the median estimates for the eclipse center, eccentricity, 
flux ratio and eclipse depth from the MCMC analysis. The un- 
certainties correspond to the 68.2% credible intervals, or the lcr 
confidence limits if the posterior distributions are close to nor- 
mal. Finally, the last two columns list the 95% credible inter- 
vals for the brightness temperature (the temperature estimated 



from the eclipse depth assuming the planet to radiate as a black 
body), 7br, and equilibrium temperature (the temperature the 
planet would have under thermal equilibrium if the stellar ra- 
diation would be the only heat source), r eq . 

The panels in the Figs. [TJ [2] [3] |4] and [5] plot (a) the model 
probability and (b) Bayes factor as a function of the eclipse cen- 
ter time; (c) the CoRoT light curve and the model derived from 
the MCMC analysis, both phase-folded and binned using three 
bin widths; and (d) the marginal posterior distributions from the 
MCMC analysis for the flux ratio, (e) eclipse center phase, (f) 
orbit eccentricity, and (g) argument of the periastron. 

We list the secondary eclipse ephemerides in Table [6] and 
show the expected flux ratio, eclipse depth and equilibrium tem- 
perature limits for planets without eclipse candidates in Table [7] 
The equilibrium temperatures are calculated assuming the planet 
to radiate as a black body, and the estimates are calculated using 
the published values for the effective stellar temperature, planet 
radius ratio and semi-major axis. 

The brightness and equilibrium temperature estimates are 
calculated in Monte Carlo fashion combining the published es- 
timates for the semi-major axis, radius ratio, and stellar temper- 
ature with our flux-ratio estimate. The equilibrium temperatures 
are calculated as 

r eq = r*0- 1/2 (/(i-A B )) 1/4 , (13) 

where is the estimated stellar temperature, 9. d the scaled semi- 
major axis, / the heat redistribution factor, and Ab the Bond 
albedo. The brightness temperature, T\, r , is numerically solved 
from the equation 

e { =A g 6 2 A +B(T br )/B(T ir ), (14) 

where the first term represents the reflected light and the second 
the thermal radiation with 6f being the surface flux ratio, A g the 
geometric albedo (we use A g = 1.5Ab), and B is Planck's law 
(i.e., we approximate the stellar radiation with a blackbody). The 
planet's Bond albedo is allowed to range from to 0.3, and the 
heat redistribution factor from 1/4 to 2/3. Thus, the reported 95% 
credible intervals advises us on the possible values of 7b, and T eq 
assuming a small albedo and without assumptions about the heat 
redistribution efficiency. 
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6. Discussion 

6. 1 . Non detections 

The Bayes factor maxima in support of the model M\ for the 
planets listed in Table [7] are all below unity, or the peaks do not 
pass the plausibility tests. In most cases the planets without de- 
tected eclipse candidates have a small expected T eq , their host 
star shows strong small-scale variability, or the host star is rela- 
tively faint. 

The only major exception is CoRoT-7b with a high equilib- 
rium temperature and a bright host star. However, the small size 
of CoRoT-7b makes the detection of its eclipse highly unlikely. 

CoRoT-14b has the highest expected equilibrium tempera- 
ture of the set, but its host star is very faint (V=16), and the 
light curve is dominated by correlated small time-scale noise. 
The Bayes factor mapping of CoRoT-14b light curve shows sev- 
eral peaks, but they can all be attributed to single events due to 
correlated noise. 




-2-1012 -20 



Time [h - 0.5 x period] Time [h - 0.5 x period] 



6.2. Previously reported eclipses 
6.2.1. CoRoT-1b 



We use CoRoT-lb (Bar ge et al.|2008) > as a first test case for our 
method. The light curve has been shown to contain an eclipse 
( [Alonso et aLp 009b; S nellen et al.|2009) , which has been con- 
firmed with ground-based near-infrared (NIR) observations by 
Rogers et al.| ( |200^ , piUorretll^ ( |2009) and |Zhao et aL] ( |20l2) . 



We identify the eclipse with an M\ posterior probability of 
88%. The results from the model comparison and parameter es- 
timation are presented in Tableland Fig. T] Our flux ratio esti- 
mate of 1% ± 0.4% yields an eclipse depth of (2 ± 0.8) x 10~ 4 , 
which is within the 1 cr limits of the white l ight curve depth of 
(1.6 ± 0.6) x 10~ 4 reported byl Alonso et al. l and the r ed channel 
depth of (1.3 ± 0.4) x 10" 4 obtained by [Snellen et al 



Compared to the monochromatic depth by Alonso et al. the 
marginally deeper transit in the red channel is expected if the 
planet has a small albedo in visible light. Our results diverge 
more strongly from those by Snellen et al.[ who find a smaller 



depth in the red channel than Alonso et a 



in the white light. 
However, the differences are not statistically significant, and can 
be due to the differences in the data reduction. 



6.2.2. CoRoT-2b 



CoRoT-2b ( |Alonso et al.||2008"] ) is used as a second test case. 
The light curve contains an eclipse signal ( |Alonso et al.|20 09b ) 
that has again later been confirmed by ground-based NIR obser- 
vations ( |Alonso et al.||2010l Caceres, C, et al., based on VLT 
observations, personal communication). 

We identify the eclipse with a marginal Mi probability of 
66% (Bayes factor of 2). Alonso et al. ( 2009b) report a secondary 
eclipse depth of (6.6 + 2) X 10" 3 at orbital phase 0.494 ± 0.006, 
values that agree with our estimates of a depth of (7 + 5) x 10~ 5 
at phase 0.501 ± 0.002. Our uncertainty estimate for the depth is 
considerably greater than the estimate by | Alonso et al. which is 
probably due to the different approaches used for the parameter 
estimation. 




_L_L 



0.5 1.0 1.5 
Flux ratio [%] 

1 



2.0 0.488 0.492 



0.496 
Phase 



0.500 




0.08 0.12 
Eccentricity 



Fig. 1. Main results for CoRoT-lb. Shown are the eclipse model 
probability and the Bayes factor (a and b, respectively); the 
lightcurve and the model, both folded over the period and binned 
with three different bin widths (c); and posterior distributions for 
the flux ratio, eclipse center phase, eccentricity and the argument 
of the periastron (d, e, f, g) 



6.3. Statistically significant eclipse candidates 
6.3.1. CoRoT-6b 

CoRoT-6b is a hot Jupiter with M P = 2.96M Jup and R P - 
1 . 177?j up orbiting an F9V star on a relatively long-period or- 



bit of 8.9 days ( |Fridlund et al^|2010| ). The light curve fea- 



tures an eclipse candidate with a small but non-zero eccentricity 
(6> e = 0.06 ± 0.01). We obtain a Bayes factor of 10.6 (see Fig. [3), 
corresponding to posterior M\ probability of 91.4%. The candi- 
date passes all our tests, and does not seem to be due to a single 
event. 
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Fig. 2. As in Fig. [I] for CoRoT-2b. 



Fig. 3. As in Fig.[l] for CoRoT-6b. 



However, the candidate is deep (f3f = 1.96% +0.69%) consid- 
ering the planet's period, and yields a minimum planetary bright- 
ness temperature that is nearly 1000 K higher than the maximum 
equilibrium temperature. If the eclipse is real, substantial amount 
of additional heating would be needed to explain the depth. 

6.3.2. CoRoT-11b 



The light curve of CoRoT-1 lb ( jGandolfi et al.|2010] > shows the 
most significant eclipse candidate identified by our study. We 
obtain a Bayes factor of (5.5 ± 1.1) x 10 4 , corresponding to a Mi 
probability of 99.998%, for an eccentric orbit with 6 e = 0.35 and 

= 15°. The candidate signal passes the tests, and is unlikely 
to be due to a single event or correlated noise. 

CoRoT-llb is an inflated planet with M P - 2.33 M Jup and 
Rp = 1.4 j?j up , orbiting on a 2.99 d period a hot rapidly rotating 
F6V star. |Gandolfi et al.| ( |2012[ ) measure a sky-projected spin- 



orbit angle of 0.1° + 2.6° using Doppler tomography, and con- 
clude that the orbit of the planet is most likely closely aligned 
with the stellar rotation axis. 

The rapid rotation of the primary makes accurate estimation 
of the eccentricity from radial velocity observations difficult. 
Gandolfi et al. ( 20 1 0) > give an RV-based constraint of f? e < 0.7, 
and comment that for 6 e <: 0.2 the mean stellar density p esti- 
mated from the transi t would be inco mpatible with an F6 dwarf 

is 0.69+0.02 g/cm 3 ). Using 



star (the p reported by Gandolfi et al., 

the equation (4) for the mean stellar densityQp given by Tingley 
[etaLl ( |2011a| l, 



30 P g J /(l +0 k f-b 2 Q 
7T 2 Gr 



14 



2^2\3/2 



1 



(15) 



7 We reproduce here a version of the equation corrected for a missing 
Q 2 term (Tingley, 2012, private communication). 
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where T14 is the transit duration, G the gravitational constant, 
and Q = (1 - + 9 P , sin ft,,), with ou r e and 9 W values, we 



obtain p = 0.51 g/cm . Gandolfi et al. state that p would be 



unrealistically high for high eccentricities, but this applies only 
for 9 W = 0. However, if we decrease the impact parameter to b = 
0.39 (necessary if the eclipse is considered to be real, see below), 
we can reproduce the reported value of p — 0.69 + 0.02 g/cm 3 , 
and p is within 3rx from its reported value for 0.15 < b < 0.55. 

The value of 9 W = 75° implies that the eclipse occurs near 
the apoastron, and the transit near the periastron. The impact 
parameter for the transit must therefore be smaller then the im- 
pact parameter for the eclipse. The current transit-derived impact 
parameter estimates of 0.8 ( Gandolfi et al.|[20Tm |Southworth 
201 1 Gandolfi et al.|2012 1 have been calculated assuming a cir- 
cular orbit, and only a dedicated light curve modeling beyond the 
scope of this work could show if the transit could be reproduced 
with the given 9 e , 9 W and a close-to-zero impact parameter. 

We derive a planet-to-star flux ratio of 3.2%, leading to 
an eclipse depth of (3.6 + 0.7) x 10~ 4 . This corresponds to 
7br = 2800 K, which is significantly higher than the maximum 
equilibrium temperature of 2200 K. The high temperature might 
be explained by tidal heating due to high eccentricity ( |Jackson| 
et al.|2008 Leconte et al.|2010 1, but the analysis of the physical 
plausibilities of different heating mechanisms is again outside 
the scope of this work. Further RV measurements to better con- 
strain the eccentricity are subject to ongoing observations. 

Given the deep eclipse in the visible light and the brightness 
of the primary (V=12.8), the eclipse should be confirmable using 
ground-based NIR observations. Rough predictions for the J, H, 
and K band depths are 1.0, 1.6, and 2.2 mmags. 

6.3.3. CoRoT-15b 



CoRoT-15b (Bouch y" et al.|2010| ) is an inflated high-mass brown 
dwarf with M = 63M Jup and R = 1 . li? Jup on a P = 3 day orbit 
around an F7V star. The star is faint (V « 16), and the light curve 
is monochromatic, contains only long-cadence data, and shows 
several jumps and strong stellar variability. 

We identify an eclipse candidate with a Bayes factor of 27, 
corresponding to a Mi probability of 97%. The candidate signal 
is deep, yielding a minimum brightness temperature estimate of 
3500 K. The temperature is suspiciously high, but given the high 
mass, inflated radius, stellar insolation, and possibly young age 
of the system (Bouchy et al. give two age estimates derived using 
different methods, 1.14—3.35 Gyr and 1.9 + 1.7 Gyr, where the 
latter method gave possible pre-main sequence solutions), may 
still be physically feasible. 



6.4. Marginal detections 

In addition to the statistically significant candidates, we also 
report four marginal candidates that pass the plausibility tests. 
Considering their Bayes factors, all the marginal candidates fall 
into the "not worth more than a bare mention" class (Jeffreys 
T96Tj|Kass & Raftery|1995[ >. 



CoRoT-3b is an object near the planet/brown dwarf boundary 
(M P = 21.66 Mj up , R P = 1.01 flj U p) in a 4.26 day or- 



bit (Deleuil et al. [2008). We find a tentative eclipse near 
an orbital phase of 0.51. The candidate is marginal with a 
Bayes factor Bio of 2.3, and the light curve shows corre- 
lated noise on the time-scale of hours. The 95% credible in- 
tervals for 7b r and T eq overlap, but the median brightness 
temperature (2700 K) is high compared to the upper limit 
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Fig. 4. As in Fig.[l] for CoRoT-llb. 



equilibrium temperature of 2200 K. However, with its large 
mass, CoRoT-3b is above the deuterium-burning mass limit 
(Spie gel et al.||20l"T >, and the additional heating from deu- 
terium burning early during the planet's history might be 
able to explain the temperature. 
CoRoT-13b is a high-density hot Jupiter with Mp = 1.3 M Jup , 
R P = 0.9 R Sup , and P = 4.04 days ( |Cabrera et al.|2010| . The 
light curve features several jumps, and the presence of low- 
amplitude small-time-scale variability is evident. We find a 
tentative eclipse passing our tests with Bio = 3.3 near a phase 
of 0.48. Unlike CoRoT-3b, the T br and r eq 95% credible in- 
tervals do not overlap, and the median (2800 K) is signif- 
icantly higher than the upper limit equilibrium temperature 
(1700 K). 

CoRoT-18b is a massive hot Jupiter with M P = 3.47 M Jllp and 
Rp = 1.31 j?j up on a short-peri od orbit (P — 1.9 d) around a 
G9V star (Heb rard et al.|20lT] >. The primary is likely to be 
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Fig. 5. As in Fig.[l] for CoRoT-15b. 



a young star with an age of several hundred Ma, but the age 
is not well constrained. We find a secondary eclipse candi- 
date near a phase of 0.47, corresponding to an eccentric orbit 
with 6> e = 0.1 ±0.04. As with CoRoT-3b and CoRoT-13b, we 
obtain a high flux ratio corresponding to a median bright- 
ness temperature (2700 K). If we were to assume the eclipse 
to be real, the high eccentricity and temperature would both 
support the young stellar age hypothesis. 
CoRoT-21b is a hot Jupiter with M P = 2.26 M Jup , R P = 
1.30 Rj up on a P = 2.72 day orbit around a faint (V=16) 
F8IV star (Pat zold et al.|2012] l. We find a marginal eclipse 
candidate with B\q = 1.7. The signal passes the splitting 
test, but the mean brightness temperature estimate of 3100 K 
is again anomalously high. However, the flux ratio estimate 
has a large uncertainty, and the 95% credible intervals of T^ r 
and r eq overlap. 



6.5. Shortcomings and possibilities for future improvement 

While our brightness temperature estimates for CoRoT-lb and 
CoRoT-2b agree with the published values, all the new candi- 
dates show anomalously high brightness temperatures. For some 
of the candidates the temperatures could be explained by young 
age, atmospheric temperature inversion, deuterium burning, or 
tidal heating. However, the plausibility of different additional 
heat sources would need detailed modeling of the individual sys- 
tems that is out of the scope of the current generalizing study 
aimed at presenting the new candidates. 

Assuming normally distributed i.i.d. errors is a simplifica- 
tion, especially since the light curves are known to contain 
sources of systematic noise. The assumption of normality was 
tested to be a good approximation for the analyzed light curves, 
but the errors cannot be considered independent. Inclusion of 
correlated noise may increase our ability to distinguish true 
eclipse signals from the noise due to the stellar granulation, and 
will be considered for the future studies. 



7. Conclusions 

We have presented and characterized three new statistically sig- 
nificant eclipse candidates and four new marginal candidates, 
and offered an independent confirmation of the CoRoT-lb and 
CoRoT-2b eclipses. 

All three of the new significant candidates have higher Bayes 
factors than the two confirmed eclipses, but they either imply 
properties that deviate from the ones derived from transits (e.g., 
the impact parameter of CoRoT-llb), or yield high brightness 
temperatures. Three of the four marginal candidates have higher 
Bayes factors than the CoRoT-2b eclipse. These candidates cor- 
respond to planets orbiting faint host stars and the light curves 
contain significant amounts of correlated noise, faintness and 
noise both contributing to the large uncertainty in the derived 
parameters. 

The estimated optical bandpass brightness temperatures are 
higher than the expected equilibrium temperatures for all of our 



new candidates. This agrees with the results by |Coughlin & 
Lopez-Morales] ( |2012| l, who find significant amounts of excess 
flux for their Kepler eclipses, and can be explained simply by the 
fact that planets are not black bodies. However, further observa- 
tions of eclipse events in the IR are necessary to confirm both 
the candidates and the characterization presented in this work. 
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